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Abstract In this paper we propose a new model for quasi-periodic oscillations (QPOs) based on oscillating 
magnetohydrodynamic modes in neutron star magnetospheres. We argue that the interaction of the accretion 
disk with the magnetosphere can excite resonant shear Alfven waves in a region of enhanced density gradients, 
the region where accretion material flows along the magnetic held lines in the magnetosphere (see Fig. 1). We 
demonstrate that depending on the distance of this region from the star and the magnetic held strength, the 
frequency of the held line resonance can range from several Hz (weaker held, farther from star), to approximately 
kHz frequencies (stronger held, ~ 2 — 10 stellar radii from the star). We show that such oscillations are able 
to signihcantly modulate inhow of matter from the high density region toward the star’s surface, and possibly 
produce the observed X-ray spectrum. In addition, we show that the observed 2 : 3 frequency ratio of QPOs is a 
natural result of our model. 
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1. Introduction 

The discovery of quasi-periodic oscillations (QPOs) in low- 
mass X-ray binaries (LMXBs) has been reported and dis¬ 
cussed in recent studies in theoretical and observational 
astrophysics. The Rossi X-Ray Timing Explorer has ob¬ 
served oscillations in the X-ray flux of about 20 accreting 
neutron stars. These oscillations are very strong and re¬ 
markably coherent with frequencies ranging from ~ 10 Hz 
to ~ 1200 Hz. The kHz frequencies correspond to dynami¬ 
cal time scales at radii of a few tens of kilometers, and are 
possibly closely related to the Keplerian orbital frequency 
of matter at the inner disk. Almost all sources have also 
shown twin spectral peak QPOs in the 500 Hz to 1200 Hz 
part of the X-ray spectrum, and both peaks are moving 
up and down in frequency simultaneously. 

Although there are some detailed differences in differ¬ 
ent types of X-ray sources, the observed QPOs are re¬ 
markably similar, both in frequency and peak separation. 
(Six of the 20 known sources are originally identihed as 
“Z” sources and the rest are known as “atoll” sources 
(see table 1). For more information about the atoll and 
Z sources, see van der Klis 2000). Such a similarity shows 
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that QPOs should depend on general characteristics of the 
X-ray sources which are common to all systems. In other 
words, the QPO can be regarded as a generic feature of 
the accreting neutron star. 

Remarkably, the characteristic dynamical time 
scale for material orbiting near the compact ob¬ 
ject is comparable to the observed millisecond X- 
ray variability, ie. because r^yn = y/r^/{GM) ~ 

2 ms (r/100km)^/^(1.4MQ/M)^/^ where r is the distance 
of the orbiting material from the center of the compact 
object with mass M (circular orbit), G — 6.67 x 10“® 
dyn cm^/g^, and Mq = 1.99 x 10^^ g. Such a natural 
time scale is the foundation of most models for the 
observed QPOs. Alpar & Shaham (1985) first proposed 
a beat-frequency model to explain the low-frequency^ 
horizontal branch oscillation (HBO) seen in Z sources 
(see also Lamb et al. 1985). They used the Alfven 
radius, where the magnetic pressure balances the ram 
pressure of the infalling material in spherical accretion, 
ta ~ 1.5 X 10® cm {M, as the 

preferred radius. For this Alfven radius, the assumptions 

^ Two different low-frequency (< 100 Hz) QPOs were known 
in the Z sources, the 6 — 20 Hz so-called normal and flaring- 
branch oscillation (NBO; Middleditch & Priedhorsky 1986) 
and the 15—60 Hz so-called horizontal branch oscillation (HBO; 
van der Klis et al. 1985). 
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are: magnetic field is stellar dipole, velocity and mass 
density are for free-fall. Here Mn is the mass accretion 
rate (this corresponds to 1.5 x lO“®M 0 yr“^ and is typical 
for low-mass accreting neutron stars) in units of 10 ^^ g 
s“^ and /i 26 is the magnetic dipole moment at the stellar 
surface in units of 10^® G cm^. Note that we assumed M 
is constant in time and position. 

Miller et al. (1996, 1998) proposed a beat-frequency 
model to study kHz QPOs, the so-called the sonic-point 
beat-frequency model, based on a new preferred radius, 
the sonic radius. In this model, they assumed that near 
the neutron star there is a very narrow region of the disk 
in which the radial inflow velocity increases rapidly as the 
radius decreases. Such a sharp transition in the radial ve¬ 
locity of plasma flow from subsonic to supersonic happens 
at the “sonic point” radius, Tsonic. This radius tends to 
be close to the innermost stable circular orbit (ISCO) 
’’ISCO) however, radiative stresses may change its location, 
as required by the observation that the kHz QPO frequen¬ 
cies vary with time. Comparing the HBO and kHz QPO 
frequencies, clearly rsonic ‘C ta, so part of the accreting 
matter must remain in near-Keplerian orbits well within 
ta. An important prediction of the sonic-point model is 
that Ajz = — be constant at the spin frequency of the 

neutron star Vg, which is contrary to observations (Lamb 
& Miller, 2001). Accurate observations from the millisec¬ 
ond pulsar SAX J 1808.4-3658 show that the frequency sep¬ 
aration = 195 Hz, while the pulsar spin frequency is 
Vg = 401 Hz (Wijnands et ah, 2003; Chakrabarty et ah, 
2003). Nonetheless, further observations reveal that the 
frequency separation A^ decreases when increases. 

Stella & Vietri (1998, 1999) proposed a relativistic pre¬ 
cession model in which the high kHz QPO frequency V 2 is 
identified with the Keplerian frequency of an orbit in the 
disk (similar to the beat-frequency model) and the low 
kHz QPO frequency vi with the periastron precession of 
that orbit. Furthermore, in the relativistic model A^ and 
Ug are not expected to be equal as in a beat-frequency in¬ 
terpretations. An interesting result of the relativistic pre¬ 
cession model is that /S.v should decrease not only when 
V 2 increases (as observed) but also when it sufficiently de¬ 
creases. However, in this model it is not clear how the 
X-ray flux is modulated at the these frequencies. 

Recently Zhang (2004) proposed a model for kHz 
QPOs based on magnetohydrodynamic (MHD) Alfven os¬ 
cillations in the disk. He introduced a new preferred ra¬ 
dius, the quasi-sonic point radius, where the Alfven ve¬ 
locity matches the orbital Keplerian velocity. The author 
suggested that the upper and lower kHz frequencies are 
the MHD Alfven wave frequencies at the quasi-sonic point 
radius with low and high accreted material mass densities, 
respectively. In this model, the frequency separation will 

^ In general relativity, no stable orbital motion is possible 
within the innermost stable circular orbit (ISCO), nsco = 
GGMjc^ « 12.5Mi, aMq km. The frequency of orbital motion 
at the ISCO, the highest possible stable orbital frequency, is 
J^isco ~ (1580/Ml,4M0 ) Hz. 


decrease when vi increases and/or sufficiently decreases. 
Although his results show good agreement with observa¬ 
tions, several vital issues are left with no answers. For 
example, no mechanism is suggested to explain the exci¬ 
tation of Alfven wave oscillations within the star’s mag¬ 
netosphere. More importantly, it is not clear how these 
Alfven wave frequencies modulate the X-ray flux coming 
from the surface of the star. 

In this paper, we propose a new model to explain QPOs 
in LMXBs based on interaction of accreting plasma with 
the neutron star magnetosphere. This interaction is a com¬ 
mon feature in all accreting binary systems. Matter trans¬ 
fers from the companion star to the compact object. This 
matter is accelerated by the gravitational pull of the com¬ 
pact object and hits the magnetosphere of the star with 
a sonic/supersonic speed. The MHD interaction of the in¬ 
falling plasma with the neutron star magnetosphere, will 
alter not only the plasma flow toward the surface of the 
star, as assumed by current QPO models, but also the 
structure of the star’s magnetosphere. The magnetic Held 
of the neutron star is distorted inward by the infalling 
plasma of the Keplerian accretion flow. Since the grav¬ 
ity of the star confines the inward flow to a small solid 
angle ~ 10 “^ ster, the magnetic field of the star will be 
more compressed in the disk plane than in other areas 
(see Fig. 1). Furthermore, in a more realistic picture, one 
would expect that the highly accelerated plasma due to 
the infall process would be able to penetrate through the 
magnetic field lines. Such material will be trapped by the 
magnetic field lines and produce enhanced density gradi¬ 
ent regions within the magnetosphere. ^ However, due to 
the material’s initial velocities, the penetrating material 
moves along the field lines, finally hitting the star’s sur¬ 
face at the magnetic poles and producing the observed 
X-ray fluxes (Fig. I). See Ghosh et al. (1977) for more 
detail. 

Besides modifying the geometry of the magnetosphere, 
the compressional action of the accretion flow can excite 
some plasma perturbations in the region of enhanced den¬ 
sity gradients. This can be understood by noting the fact 
that the inward motion of the accretion flow will be halted 
by the outward magnetic pressure at a certain distance 
from the star, the Alfven radius (Ghosh et ah, 1977). In 
other words, the accretion flow pushes the stellar mag¬ 
netic field toward the star until the pressure of the infalling 

® Note that the enhanced density gradient region is defined 
as follow: when we move outward from the surface of the star, 
the plasma density in the magnetosphere is more or less is con¬ 
stant, say pQ. At the point that magnetosphere interacts with 
the accretion disk (more likely at the disk’s inner edge), how¬ 
ever, there is a sudden increase in the magnetospeheric plasma 
density, say p > po, due to the motion of accretion materials 
along field lines. After this thin region (comparing with whole 
magnetosphere), the magnetospheric plasma density falls down 
to po again. Therefore, the enhanced density gradient region is 
not a region with a gradient in the plasma density, only. It is 
a region with different (higher) plasma density compared with 
the rest of the magnetosphere. 
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plasma pv "^/2 balances the magnetic pressure Here 

p and Vj. are the density and radial velocity of the in¬ 
falling matter and Bp is the poloidal magnetic field at 
the disk plane. Therefore, any instability in the disk at 
the Alfven radius would disturb such quasi-equilibrium 
configurations as well as the structure of the magneto¬ 
sphere, possibly through the Rayleigh-Taylor and/or bal¬ 
looning instability. As an example, the interaction of the 
solar wind with the Earth’s magnetosphere excites reso¬ 
nant shear Alfven waves, or field line resonances (FLRs), 
along the magnetic field lines (Samson, 1991). As a result, 
one might expect that such Alfven waves can be excited by 
the compressive accreting plasma in the magnetosphere of 
an accreting neutron star with the accretion flow playing 
the role of the solar wind. 

The plan of this paper is as follow. We outline the oc¬ 
currence of resonant coupling in a magnetized plasma in 
section 2. We review the excitation of shear Alfven waves 
by studying linear perturbations in MHD. Next, the occur¬ 
rence of FLRs, as a result of a resonant coupling between 
the compressional and shear Alfven waves, is discussed. 
To illustrate the basic features of FLRs, the excitation of 
these resonances in a rectilinear magnetic field is consid¬ 
ered. In section 3, in order to consider a more realistic 
model for an accreting neutron star, we study the influ¬ 
ence of the ambient flow along the field line on the ex¬ 
citation of FLRs. We demonstrate that in this case the 
eigenfrequency of the Alfven modes is modulated by the 
velocity of the field aligned plasma flow. Depending on the 
field line where the resonance occurs, the eigenfrequency 
of the FLR is in the range of several hundred Hz to kHz. 
Furthermore, in the presence of this flow, the plasma dis¬ 
placement parallel to the magnetic field lines is non-zero. 
This is very important because these displacements, which 
are missing in the case of zero ambient flow, might be re¬ 
sponsible for modulating the motion of infalling material 
toward the magnetic poles and producing the observed X- 
ray fluxes. Existence of this non-zero parallel component 
distinguishes our study from other related investigations 
in which the modulation of X-ray flux in observed QPOs is 
fairly unclear or ambiguous. A possible occurrence of more 
than one peak in the power spectrum is also discussed. In 
addition, we show that the observed 2 : 3 frequency ratio 
of QPOs is a natural result of our model. Section 4 is de¬ 
voted to summarizing our results and further discussions. 


2. The magnetospheric model 

2.1. Magneto-hydrodynamic waves 

In general, the dynamics of an ideal magnetized plasma is 
described by plasma density p, plasma pressure p, gravi¬ 
tational potential d>, velocity vector v and magnetic field 


B: 


dv/dt V • Vv = —Vp/p — -— B X (V X B) — V<I>, 

47rp 


dp/dt -I- V • (pv) = 0, 
dB/dt = V X (v X B), 

V-B = 0. 


(la) 

(lb) 

(l c) 

(l d) 


An adiabatic equation of state with an adiabatic index 7 , 
p/p'^ =const., is assumed in this paper in order to com¬ 
plete the above equations. 

In a steady state (ie. d/dt = 0), Eqs. (1) have been 
studied in detail in connection with the problem of stellar 
winds from rotating magnetic stars (Mestel, 1968a,b) and 
in connection with diffusing/flowing plasma into magneto¬ 
spheres in accreting neutron stars (Eisner & Lamb, 1977, 
1984; Ghosh et ah, 1977). Obviously, in the case of ac¬ 
creting neutron stars, v represents the inflow velocity of 
matter accreted to stars, while in stars with a stellar wind 
it represents the plasma outflow velocity. 

Consider an axisymmetric system consisting of a ro¬ 
tating star with a constant angular velocity along the 
z-axis and a magnetic field whose symmetry axis is aligned 
with the rotation axis of the star. By decomposing the ve¬ 
locity and magnetic field vectors into poloidal and toroidal 
components: 

V = Vp -I- B = Bp -I- B^^, (2) 


one can obtain 


Vp = (//p)Bp, (3) 

n = ns-\-if/p){B^/w), (4) 

where / is the mass flux along a magnetic flux tube of unit 
flux (Ghosh et ah, 1977). The subscripts p and (j) denote 
poloidal and toroidal components, respectively, H is the 
angular velocity of the plasma, w is the distance from the 
axis of the rotation, and </> is a unit toroidal vector. We 
note that both / and Hg are constant along a given field 
line. Equation (4) can be rewritten as 

Vl = Q.s + {vp/w){B^/Bp), (5) 

where Vp = Vp • p is the magnitude of velocity along 
the poloidal magnetic field Bp with magnitude Bp and 
p = Bp/Bp. Furthermore, Ghosh et al. (1977) have shown 
that at distances close to the Alfven radius, the magnitude 
of the poloidal component of the inflow velocity (for ac¬ 
creting systems) Vp approaches the poloidal Alfven speed 
VA, i.e. * 

_ vU^A)=vl{rA). ( 8 ) 

^ In order to estimate the poloidal component of the inflow 
velocity Vp, one needs to integrate the momentum Eq. (la): 

(l/2)(np-|-H^ro^)—GM/r = constant along a given field line, 

( 6 ) 

where M is the mass of the neutron star (Mestel, 1968a,b; 

Ghosh et ah, 1977). In the above equation the pressure term 
is neglected. Equation (6) shows conservation of energy in a 
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The possible perturbations of a magnetized plasma 
and MHD waves are found by specifying the equilibrium 
configuration of the star and then solving the perturbed 
Eqs. (1). This is a nontrivial problem, that is addressed 
by several investigators. However, in this paper we will be 
interested in the propagation of shear Alfven waves in the 
star’s magnetosphere and the resulting FLRs. 


FLRs in dipole magnetic fields can be found in Taylor & 
Walker (1984). 

To illustrate the resonant coupling between the fast 
compressional mode and the shear Alfven mode we study 
the excitation of FLRs in an inhomogeneous plasma in a 
rectilinear magnetic field model. 


2.2. Field line resonances 

In MHD, waves can propagate in three different modes 
including shear Alfven waves, and the fast and slow com¬ 
pressional modes (Landau & Lifshitz, 1992). In a homoge¬ 
neous plasma, one can easily show that these three modes 
are linearly independent. However, in an inhomogeneous 
medium these three modes can be coupled, yielding ei¬ 
ther a resonant coupling (Southwood, 1974; Hasegawa, 
1976), or an unstable ballooning mode (Ohtani & Tamao, 
1993; Liu, 1997). FLRs result from the coupling of the 
fast compressional and the shear Alfven modes ® whereas 
the ballooning instability results from the coupling of the 
slow compressional and the shear Alfven modes. However, 
in the cold-plasma approximation (which is appropriate 
to neutron star magnetospheres), the slow compressional 
wave does not exist. FLRs are standing waves that are 
stimulated in a pulsar with gradients in the Alfven speed 
transverse to the ambient B-field, particularly, within den¬ 
sity boundary layers (a parallel gradient in the Alfven 
speed is also included). Ffficient coupling between the 
shear Alfven wave and the fast compressional wave can 
produce a relatively narrow FLR spectrum, even when 
the driver is broadband. 

The linear theory of the FLRs was developed by Chen 
& Hasegawa (1974) and Southwood (1974), and applied 
to auroral phenomena by Hasegawa (1976). Samson et al. 
(2003) developed a nonlinear model with a nonlocal elec¬ 
tron conductivity to explain the evolution of field aligned 
potential drops and electron acceleration to form auroral 
arcs. They studied the possible coupling between the fast 
compressional mode and the shear Alfven mode in an in¬ 
homogeneous plasma with a radial gradient in the Alfven 
velocity va = B For further discussion of mode 
conversion for shear Alfven waves see for instance Stix 
(1992), and for an example of computational models see 
Rickard & Wright (1994). Analytic solutions for toroidal 


frame corotating with the star, while in a nonrotating frame 
the extra term appears, that represents the work done 

by the magnetic field on the flowing plasma. However, as ar¬ 
gued by Ghosh et al. (1977) the poloidal velocity in the inner 
magnetosphere where r <C ca is nearly equal to the free-fall 
velocity, i.e. 

Up ~ (2GM/r)^/^ (7) 


® In the cold-plasma approximation, the fast compressional 
wave (fast magnetosonic wave) is also called the compressional 
Alfven wave. 


2.2.1. FLRs in a rectilinear magnetic field 

Plasma dynamics can often be described by the adiabatic 
ideal MHD equations: 


p(5/at-hv-V)v = -Vp-—Bx (VxB), (9a) 

47r 

aB/5t = V X (v X B), (9b) 

V • B = 0, (9c) 

dp/dt -I- V • (pv) = 0, (9d) 

d{pp~^)/dt = Q. (9e) 


where we neglect the effect of gravitational attraction on 
the plasma (see Eqs. (1) for comparison). Introducing 
Eulerian perturbations to the ambient quantities in the 
form 


p = Po + 5p, p = po + dp, 

B = Bo-f<fB, V = Sv = d^/dt, 


we can derive a linear wave equation by simplifying Eqs. 
(9) up to the hrst order in the perturbations. We choose 
coordinates such that the ambient magnetic Held is in the 
z-direction. Without loss of generality, we assume that the 
gradients in all the unperturbed quantities are in the x- 
direction. Setting the plasma displacement as 


=^{x)e-^Gt-kyv-k.z)^ 

(10a) 

dB(r,t) 

(10b) 

Sp{r,t) = 

(10c) 

Spir,t) = 

(lOd) 


Eqs. (9) reduce to (Harrold & Samson, 1992) 

, F'{x)d^, 


dx'^ F{x) dx 


G{x)^x = 0 , 


with 

G(x) = 
and 


( 11 ) 


^ ^ - (vl + vl){kl + fc2) + kl{kl + ki)vlvl 


F{x) = 


{vl + vDuj"^ - k^vlvl 
Po(w^ - klv\) 


G{x) 


( 12 ) 

(13) 


where v\{x) = i?Q/(47rpo) and Vg{x) = jpojpo, and 
' = d/dx. Here we assumed that the total pressure 
(fluid-l-field) Po = Po + Bq/{8tt) is constant. Equation (11) 
yields two turning points at G{x) = 0, and two resonances 
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at F{x) = 0. Close to the resonance positions, Eq. (11) can 3. FLRs in the presence of ambient flow (along 
be approximated by the magnetospheric field lines) 




1 


X — xo dx 


+ G{x)^x — 0 , 


(14) 


where xq is the position where the resonance occurs. 

In the case of a strong magnetic field, or for a cold 
plasma (ie. Pfluid/Pmagnetic < 1), one can put us - 0. 
Then G(x) and F{x), Eqs. (12) and (13), reduce to 

G(x) = '^-kl-kl (15) 


and 


^ - klv\) 

- (k^ + k‘l)v\ ' 


(16) 


In this case Eq. (11) has only one resonance (the Alfven 
resonance) at 

- klvl = 0, (17) 


which corresponds to the dispersion relation for the shear 
Alfven wave along the field line. It is necessary to note 
that in above calculations four boundary conditions are 
assumed. In the x-direction the boundaries, say, at Xmin = 
0 and Xmax = Xm are considered perfectly reflecting (i.e. 

= 0 at boundaries). Eurther, we assume that neutron 
star and accretion disk surfaces are perfectly conducting 
so that the displacement ^ vanishes at these surfaces. In 
the rectilinear model we represented these boundaries by 
planes at constant z, i.e. Zmin > 0 on the disk and Zmax > 
0 on the star’s surface (see Fig. I). 

Briefly, the FLR mechanism can be outlined as follows: 
the incoming compressional wave hits the stellar magneto¬ 
sphere and reaches a field line (or a magnetospheric shell) 
in which the frequency of the incoming wave matches the 
eigenfrequency of the standing shear Alfven wave along 
that particular filed line. The resulting forced-oscillating 
system causes the amplitude of the shear Alfven wave to 
grow in time forming FLRs. 

The FLR mechanism is generic and likely to occur 
in many astrophysical magnetospheres. As a result, one 
would expect that FLRs likely occur not only in the 
Earth’s magnetosphere but also in the magnetospheres of 
accreting neutron stars. In the case of the Earth’s magne¬ 
tosphere, the source of energy for the FLRs is the inter¬ 
action of the solar wind with the magnetosphere (Harrold 
& Samson, 1992). In accreting neutron stars, the accreted 
plasma interacts with the stars’ magnetosphere, allowing 
the compressional mode to propagate into the magneto¬ 
sphere and flow along the field lines toward the magnetic 
poles. Such a compressional action of the accretion flow 
can excite resonant shear Alfven waves in the enhanced 
density regions filled by plasma flowing along the field 
lines. In section 3 we consider this mechanism more care¬ 
fully to address its possible relation to the QPOs observed 
in accreting neutron stars in LMXBs. 


As discussed above, FLRs have been used to model elec¬ 
tron acceleration and auroral arcs in the Earth’s magne¬ 
tosphere. Although FLRs are likely to be excited in any 
magnetospheric system with an input of compressional en¬ 
ergy, one must carefully evaluate the differences between 
the Earth’s and neutron star magnetospheres. In neutron 
star magnetosphere, for example, the strong magnetic field 
of the neutron star, the rapid stellar rotation, and the in¬ 
tense radiation pressure from the stellar surface should be 
considered. 

In the case of the Earth’s magnetosphere, due to the 
small gravitational attraction of the Earth and also its 
large distance from the Sun, the solar wind impacts the 
whole Sunward side of the geomagnetosphere. This super¬ 
sonic solar wind produces the so-called bow shock struc¬ 
ture at the outer boundary of the geomagnetospher. On 
the contrary, the strong gravity of the neutron star cre¬ 
ates a converging flow that is supersonic long before the 
flow hits the star’s magnetosphere. Such a localized flow 
is able to change the structure of the magnetosphere in lo¬ 
cal areas, particularly in the equatorial plane. In addition, 
the highly variable nature of the exterior flow can change 
the magnetosphere’s structure dramatically in time. Large 
flux of plasma stresses the star’s outer magnetosphere and 
creates a relatively high plasma density in this region. The 
plasma then flows along the field lines, an interior flow, un¬ 
til it hits the star’s surface near the magnetic poles (see 
Fig. 1). In this section we study the excitation of FLRs 
by considering such plasma flows in the neutron star mag¬ 
netosphere. In the context of the Earth’s magnetosphere, 
however, Harrold & Samson (1992) studied the resonant 
excitation of Alfven waves by surface waves excited by the 
solar wind and plasma flows in the Earth’s magnetosheath. 

The presence of a flow v in the plasma adds more 
modes to the plasma waves. In general, such flows are a 
combination of plasma flow along the magnetic field lines, 
Vp, and rotational motion of the plasma around the star, 
with angular velocity ft, ie. v = Vp + Clx r (see Eq. 2). As 
we shall see in Eq. (22a), the existence of Vp alone induces 
a non-zero parallel component (relative to the direction of 
the magnetic field) of displacement that vanishes in a cold 
plasma if Vp = 0. We note that this component can be 
responsible for modulating the infalling plasma flow near 
the star’s surface with FLRs’ frequencies and producing 
the observed X-ray fluxes, see Eq. (22a) below. We will 
return to this point later. 

In a strong magnetic field, the usual definition of the 
Alfven velocity, v\ = /{Attp), is not correct as with this 
non-relativistic definition, velocities can be larger than the 
speed of light, c. The correct definition, the relativistic 
Alfven velocity, uses energy density e relative to c^, which 
for an ideal gas reads 

e _ -I- 7 p /(7 — I) -I-R^/(47r) 


(18) 



6 


Rezania & Samson: Resonant Alfven waves in neutron star magnetospheres 


rather than the plasma density p. Consequently the (rel¬ 
ativistic) Alfven velocity is 


(19) 

In Fig. 2, we plot radial profiles of the relativistic Alfven 
velocity va in the equatorial plane (assuming p = 0). As 
shown in this figure, close to the star the relativistic Alfven 
velocity approaches the speed of light, ie. ua ~ c. This def¬ 
inition for the Alfven velocity is used extensively in stud¬ 
ies of relativistic jet dynamics in black hole systems and 
quasars (see for an example Mizuno et al. 2004). Again, we 
use a linear approximation to obtain a dispersion relation 
with resonant coupling. 

The linearized perturbed magnetohydrodynamic equa¬ 
tions in the presence of an ambient flow can be obtained 
from Eqs. (9) as 


( ^ + V • V,5v -h (5v Vv ) = -V5P 
\ at 


VP 


6h 


( 

47r V 


dB VB + B VSB 


dSB 


= V X ((5v X B -I- V X <5B), 


dt 

V • (5B = 0, 
d 

(— -I- V • V)5p + dv • Vp = —"f(6pV ■ V + pV ■ dv), 
at 


(20a) 

(20b) 

(20c) 


(20d) 


where Sv = d^/dt. P = p + i?^/(87r) and dP = Sp + B ■ 
dB / (47r) are the total unperturbed and perturbed plasma 
pressure (fluid -|- field), respectively. Here h,p,v, and B 
are the unperturbed quantities while dh, 6p, Sv, and dB 
are the perturbed quantities. Furthermore, we consider the 
slow rotation approximation and neglect both the toroidal 
field B(j, and velocity Eq. (2), to avoid complexities. 
Such assumptions may not meet the actual configuration 
precisely. Nevertheless, these assumptions simplify our cal¬ 
culations significantly. 

To analyze the problem analytically, we consider again 
a rectilinear magnetic field configuration. We note that al¬ 
though this configuration may not be suitable for the ac¬ 
creting neutron star, it provides us with a descriptive pic¬ 
ture that can be applicable to QPOs. Separating Eqs. (20) 
into parallel and perpendicular components relative to the 
ambient magnetic field and assuming perturbed quantities 
in the form of 


(j(r,t) =(j(r_L)e-*(“*-'=iini), (21) 


we find 


-tu;D^\\ = -^V ( 22 a) 

UJd 

( 2 7,2 2 _IjJdV±6P LOd VaP,. 

{uJD-kiiV^)^j_ = — - - --— j^Sh, ( 22 b) 

—iojoSB^^ = iu!B{V± ■ ^a) + *i^(^a ' ^±)B 

-I-((IBa • Va)u|| , (UsingV • (5B = 0) (22c) 

—iuJDSB± = ujk\iB^j_, (UsingV • B = 0) (22d) 

-iuJuSp = zw(|a • Va)p - 7u;p(fc|j^ll - zVa • lA).(22e) 


h 


UJd 


-{v! + vl)V-^- 


icokw 


UJd 


UJd 


(22f) 


where ujd = ^~k\\v\\ is the Doppler shifted frequency, uj is 
the eigenfrequency, /cjj is the wave number in the parallel 
direction, Vp = v\\B/B, Vs = \fjpjh is the sound veloc¬ 
ity Note that we assumed that all ambient quantities 
are function of ta only. In order to derive Eq. (22), Eq. 
(22a) uses Eqs. (22d) and (22e), Eq. (22b) uses Eq. (22d), 
and Eq. (22f) is a combination of Eqs. (22c), (22d), and 
(22e). Equation (22a) shows that the plasma displacement 
parallel to the ambient field does not vanish even in the 
cold plasma approximation (p = 0)^. The non-zero ^|| can 
affect and then modulate the motion of plasma along the 
field lines. This modulation will occur at the frequency of 
the resonant shear Alfven waves uj and should be observed 
in X-ray fluxes. 

We set the magnetic field in the z-direction with gra¬ 
dients in the ambient parameters in the x-direction (rA) 
only, i.e. r’||(rA) = Vp{x) = Vz(x), and further 


(I(r, t) = . (23) 


Here we assumed that the perturbed quantities have sinu¬ 
soidal dependencies in both y- and z-directions. Equations 
(22a), (22b) and (22f) in the cold plasma approximation 
{p = 0) can then be combined into one differential equa¬ 
tion 

^ UJpE'jx) + k^v'pEjx) 
dx"^ ujd’^{x) dx 

k,VpE\x) ^ ujI, - (fcg -k Pjvj 

ujd‘^{x) \ ujd ) v\ 

(24) 




where ujd = uj — kzVp and 


F{x) = ujd^{x) 


^D - kivi 


UJ 


D 


-(kl 


k1)v 


h v\. 


(25) 


® In order to obtain Eq. (22f) we assumed that at equilibrium 
VP — V(p-|-B^/(87r)) = 0. This is valid for a non-rotating or 
very slowly rotating star. 

^ In case of zero ambient flow, n|| =0, the parallel displace¬ 
ment ^11 vanishes in the cold plasma 
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Note that here Vp = Vp{x) and va = va{x) are functions 
of X. For Vp = 0 Eqs. (24) and (25) reduce to the one 
obtained in section 2.2.1, Eqs. (11), (15) and (16) with pq 
replaced by h, as expected. 

We note that mathematically our result is similar to 
one obtained by Harrold & Samson (1992), however, they 
considered a different model in which flow in the Earth’s 
magnetosheath is perpendicular to the B-held and there is 
no ambient flow in the region of FLRs. For a constant flow 
along magnetic field lines, i.e. Vp{x) = 0, Eq. (24) reduces 
to 


F'{x) - {kl + kl)v\ 

dx"^ F{x) dx v\ 


(26) 


Equation (26) has a singular point at E = 0. The condition 
F = 0 yields the following modes: 


= k;,{vp±VA)- (27) 


The resulting modes are combination of a compressional 
mode kzVp due to the plasma flow and a shear Alfven 
mode kzVA- They produce two frequency peaks in the 
spectrum that vary with time if the velocities change with 
time. The value of kz ~ tt/L where L is the length of 
the Held line and is on the order of the radius of the 
neutron star, L ^ Rg. Therefore, for Rg ~ 10® cm and 
■^A ~ ~ 1-5 X 10® cm s“^ we get ~ 1000 Hz, which 

is comparable to the frequencies of the observed QPOs. 
In the limit Vp <C va these modes become a single mode 
of a regular Alfven resonance. In the case of ua 'C c this 
mode becomes consistent with a single mode (17). For a 
superalfvenic flow, ie. va 'C Up, however, the MHD Alfven 
waves and the resulting FLRs are suppressed due to propa¬ 
gation of the MHD wave, kzVp. In the case of superalfvenic 
plasma flow, one has 

hk‘^1? 

- t.2 1 ^2 ’ as w ^ kzVp , (28) 

which shows that Eq. (26) is no longer singular. Therefore, 
no FLRs are likely to occur in this case. 

As a result, we expect that the two FLR peaks (Eq. 
27) occur where Vp and va are more or less comparable 
(i.e. at the same order of magnitude). As we now show, 
using the flow velocity and the Alfven velocity definitions 
below, one might expect that the FLRs likely occur at 
Rg < r < 100 i?s- 

Approximating the plasma inflow velocity with the free 
fall velocity Vp ~ vs{r) and va ~ B{r)/\/4:TTps where 
ps = M/{ve 47rr^) is the free fall mass density, one can 
rewrite Eq. (27) as: 

oj^{r) ~ kz[^ ± B(t) j^psij) ^ , 

~ kz(vp{Rg)x~^l'^ ±VA{Fig)x~'^^^), (29) 

where Vp{Rg) = {2GM / Rg)^/'^ ~ 1.6 x 

lO^®(M/M 0 )^/^(i?s/lO km)“^/^ cm s“^ and 
va{Rs) = B{Rg)/^-iTTpsiRs) ^ 4 X 


cm s“^ are 

inflow and Alfven velocities at the surface of the star®. 
Here p 2 Q is the magnetic Held dipole moment at the 
surface of the star in units of 10^® G cm®, Miy is the 
mass of accretion rate in units of 10®”^ g s“®, x = r/Rg, 
and Rg is the radius of the star. Figures 3 and 4 show 
the regular Alfven resonance frequencies {vp = 0) or 
FLRs, as a function of distance from the star and the 
Alfven speed. It is clear that the closer to the star and/or 
the larger the Alfven velocity the higher the frequency. 
Therefore, for an Alfven speed say VA{Rg) ~ -Ic one can 
get the frequency about 1000 Hz at r ~ 2.78i?s as we 
expected. Furthermore, since the Alfven velocity depends 
on the energy density (or density) of the plasma which 
varies from time to time due to several processes such as 
magneto-turbulence at boundaries during the accretion, 
the resulting FLRs frequencies will also vary with time. 
In addition, the position where the resonance takes 
place, is also subject to change in time due to the time 
varying accretion rate. Such time varying behavior causes 
time varying frequencies or quasi-periodic frequencies as 
observed. 

Setting uj~ as the observed peak separation frequency 
Av ~ uj~, and co^ as the upper QPO frequency 1^2 ~ W*", 
one can And vi, 1^2 and Av as 

vi{x) ~ 5kzVp{Rg)p2eM^^^^‘^{M 

X (i?s/I0 km)”^^"* (31a) 

V2{x) ~ kzVp{Rg)[x~^/'^ + 

2Ap26M~^^‘^ {Rg/IG km)”'^/‘‘x-®/4](31b) 

Av{x) ~ kzVp{Rg)[x~^^^ + 

2Ap2^M~^''^ {M/MqY^!^ (i?s/10 km)(31c) 

Figures 5 and 6 show the variation of Av vs v\ and V 2 ^ 
respectively, along with the observed values. The curves 
are drawn by assuming /X 26 = -32 with (top to bottom) 
kzVp(Rg) ~ 860 Hz, 750 Hz, and 650 Hz, respectively. As 
shown by observations, the value of Av decreases when¬ 
ever the magnitude of V 2 decreases (sufficiently) and/or 
it increases. Such behavior is expected in our model (see 
Fig. 6). 

® We note that to calculate the Alfven velocity at the surface 
of the star properly, one would need to use the relativistic 
Alfven velocity as defined in Eq. (19), i.e. 

VA{r) = ^ ' (30) 

^1 -h p2M-i(2GM)i/V-9/2/c2 

Here we assumed p — 0. Using the above relation the Alfven 
velocity at the surface of the star will be ~ 0.8c, .997c, and c 
for p = 10®®, 10®^, and 10®® G cm®, respectively. At further 
distances from the star, however, the relativistic Alfven veloc¬ 
ity reduces to its classical version which we used instead. The 
values given here for the relativistic Alfven velocity at the stel¬ 
lar surface, are for a neutron star with Rg = 10 km, M = 10^^ 
g s“^ and M = Mq 



Rezania & Samson: Resonant Alfven waves in neutron star magnetospheres 


The lower QPO frequency can also be rewritten as 
vi = V 2 —^v ~ 2kzVA- As a result, we obtain the frequency 
ratio vi to 1^2 

i'i/v 2 = 2va/{vp + va) = 2/(1 + Vp/ua) ■ (32) 

It is clear that depending on the relation between the 
plasma flow velocity Vp and the Alfven velocity ua, one 
can obtain different frequency ratios. For example, if the 
velocity of the inflow plasma is 2 (7/3 ~ 2.3) times greater 
than the Alfven velocity, i.e. Vp = 2 va (vp = 7/3 va), 
one can get 2 : 3 (3 : 5) ratio. For a free falling plasma 
along the lines of a dipolar magnetic field, Vp{r) oc 
whereas B(r)/oc the condition Vp ~ 2va 

(vp ~ 2.3ua) is satisfied at r ~ 3i?s ~ 1.5rA (here ^26 = 1 
is assumed). The occurrence of FLRs is very likely to be 
at this distance from the star. 

Such pairs of high frequencies with a 2 : 3 frequency 
ratio have been observed in the X-ray flux of some neu¬ 
tron stars in LMXBs® and black hole systems (the 3 : 5 
ratio has seen in one source), and it would appear that his 
feature is common to those systems. However, recent 
statistical analyzes by Belloni et al. (2004, 2005) 
showed that although the two high and low kHz 
QPO frequencies in LMXBs are well correlated, 
the frequency-frequency correlation is significantly 
different from a 2 : 3 relation. They analyzed all 
published values of kHz QPO frequencies from 
neutron star systems and found that the strong 
linear correlation observed between the lower and 
upper kHz QPO frequencies in both Atoll and Z 
sources is not compatible with a single constant 
ratio. Furthermore, by analyses of a much larger 
RXTE/PCA dataset for Sco X-1, they showed that 
there is no sharp concentration around a 2 : 3 ratio, 
but that the ratios are broadly distributed over the 
range 820 — 1150 Hz. 

The existence of such rational ratios is still a mys¬ 
tery. In black hole systems, it has been suggested that 
such frequencies correspond to a trapped g-mode or c- 
mode of disk oscillation in the Kerr metric, see for exam¬ 
ple Kato (2001) and references therein. In neutron star 
systems they are explained as the fundamental and the 
hrst harmonic of the non-axisymmetric (m = 1) g-mode 
(Kato, 2002, 2003). Further, Rezzollaet al. (2003a,b) stud¬ 
ied small perturbations of an accretion torus orbiting close 
to the black hole and modeled the observed high QPO 
frequencies with basic p-modes of relativistic tori. They 
showed that these modes behave as sound waves trapped 
in the torus with eigenfrequencies appearing in sequence 
2:3:4:... Abramowicz et al. (2003) also proposed that the 
observed rational ratios of frequencies may be due to the 
strong gravity of the compact object and a non-linear reso¬ 
nance between radial and vertical oscillations in accretion 
disks. 

® In Sco X-1 the correlation line between two frequencies has 
a steeper slope than 2/3 (see Belloni et al. 2004). 


4. Discussion 

In the present work, we have studied the interaction of 
an accretion disk with a neutron star magnetosphere in 
the LMXBs. The recent extensive observations reveal the 
existence of quasi-periodic oscillations in the X-ray fluxes 
of such stars. These oscillations, with frequencies rang¬ 
ing from 10 Hz to 1200 Hz, have been the subject of sev¬ 
eral theoretical and observational investigations. Based on 
theoretical models for the observed aurora in the Earth’s 
magnetosphere, we have introduced a generic magneto- 
spheric model for accretion disk-neutron star systems to 
address the occurrence and the behavior of the observed 
QPOs in those systems. In order to explain those QPOs 
consistently, we consider the interaction of the accreting 
plasma with the neutron star magnetosphere. Due to the 
strong gravity of the star, a very steep and supersonic flow 
hits the magnetosphere boundary and deforms its struc¬ 
ture drastically. Such a plasma flow can readily excite dif¬ 
ferent MHD waves in the magnetosphere, including shear 
Alfven waves. 

In the Earth’s magnetosphere, occurrence of aurora is 
a result of the resonant coupling between the shear Alfven 
waves and the fast compressional waves (produced by the 
solar wind). These resonances are known as FLRs and are 
reviewed in detail in section 2 of this paper. We argue 
that such resonant coupling is likely to occur in neutron 
star magnetospheres due to the interaction with accret¬ 
ing plasma. In the context of QPOs in neutron star mag¬ 
netospheres, we formulated an improved FLR model by 
considering a plasma flow moving with velocity Vp along 
the magnetic field lines. Such flows are likely to occur in 
a neutron star magnetosphere (Ghosh et ah, 1977). For a 
simple geometry, a rectilinear magnetic field, and in the 
presence of a plasma flow we found: (a) two resonant MHD 
modes with frequencies = kz{vp ± va)- (b) the re¬ 
sulting frequencies for k^ ~ tt/Rs ~ 3 x 10“® cm“^ and 
typical flow and/or Alfven velocity ~ .Ic at the stellar 
surface will be in kHz range within a few stellar radii. 
Our results match the kHz oscillations observed in the X- 
ray fluxes in LMXBs. As shown in hgures 3 and 4, the 
closer to the star and/or the larger the Alfven velocity 
the higher the frequency, (c) the quasi-periodicity of the 
observed oscillations can be understood by noting that 
due to several processes such as magneto-turbulence at 
boundaries and the time varying accretion rate, the FLR 
frequencies may vary with time, (d) a non-zero plasma 
displacement along the magnetic field lines ^||. Such a dis¬ 
placement, which oscillates with resulting frequencies 
modulates the flow of the plasma toward the surface of the 
neutron star. As a result, the X-ray flux from the star will 
show these frequencies as well, (e) setting Au = uj~ and 
r'2 = W*", one can explain the behavior of the peak separa¬ 
tion frequency Aiy relative to the upper QPO frequency ^2- 
As observed, the value of Ai/ decreases as the magnitude 
of 1^2 decreases and/or increases. Figure 6 clearly shows 
such behavior, (f) for Vp ~ 2va (obtained at r ~ 3Rs for 
M26 ~ 1)) the frequency ratio h'i/h '2 is comparable with the 
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observed frequency ratio 2 : 3 in some neutron star sys¬ 
tems. However, a broad range of frequency ratios 
(Belloni et al. 2004, 2005) is expected as values of 
Vp and VA change. 

Interestingly, using the observed values of QPO fre¬ 
quencies, one can determine the mass density and the 
magnetic dipole moment of the star using 

kzVp{r) = ( 1 / 2 )( iz 2 -I- Aiz) = V 2 - Jzi/2, (33a) 

^zfA(?") = ( 1 / 2 )( jz 2 - A^) = izi/2 . (33b) 

Therefore, 

Us ^ (327r‘‘G)"^/'‘i?®/2M^/2(M«/i?3)-i/42.i3/4^^^ 

where kz — 7r/r, /is = BqR^,x = r/Rg, and Bq is the 
magnetic field strength at the surface of the star. In table 
I, we calculate the average density Mg/R^ and magnetic 
dipole moment fig of the star for a fixed value x = 10. Note 
that larger frequencies may occur at smaller distances, ie. 
X < 10. Our results are compatible with realistic neutron 
star parameters. 

Furthermore, our model is able to explain the low 
frequency (~ 10 Hz) quasi-periodic oscillations observed 
in the rapid burster such as MXB 1730-335 and GRO 
J1744-28 (Masetti et ah, 2000). Such low frequencies can 
be extracted from the model by considering smaller in¬ 
flow/Alfven velocities and/or further distances from the 
star. 

Nevertheless, in order to avoid a number of complex¬ 
ities in our calculations, we used approximations such as 
slow rotation and a cold plasma. These assumptions may 
put some restriction on the validity of our model and re¬ 
sults. Future studies will be devoted to overcoming those 
restrictions. 
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Table 1. Observed frequencies of kilohertz QPOs in Z and 
atoll sources 


Source 

Type 

(Hz) 

1^2 

(Hz) 

Av 

(Hz) 

MJRi 

(10^^ g cm-’*) 

( 10 ’*® G cm*) 

ScoX-l 

z 

565 

870 

307±5 

2.6 

6.6 

GX5-1 

z 

660 

890 

298±11 

2.7 

7.6 

GX17+2 

z 

780 

1080 

294±8 

3.6 

8.4 

GygX-2 

z 

660 

1005 

346±29 

3.5 

7.2 

GX 340-1-0 

z 

565 

840 

339±8 

2.6 

6.6 

GX 349-1-2 

z 

710 

980 

266±13 

2.9 

8.0 

4U 0614+09 

atoll 

825 

1160 

312±2 

4.1 

8.6 

4U 1608-52 

atoll 

865 

1090 

225±12 

3.3 

9.5 

4U 1636-53 

atoll 

950 

1190 

251±4 

3.9 

10.0 

4U 1702-43 

atoll 

770 

1085 

315±11 

3.7 

8.2 

4U 1705-44 

atoll 

775 

1075 

298±11 

3.6 

8.3 

4U 1728-34 

atoll 

875 

1160 

349±2 

4.3 

9.0 

KS 1731-260 

atoll 

900 

1160 

260±10 

3.8 

9.5 

4U 1735-44 

atoll 

900 

1150 

249±15 

3.7 

9.6 

4U 1820-30 

atoll 

795 

1075 

278±11 

3.5 

8.6 

AqlX-1 

atoll 

930 

1040 

241±9 

3.1 

10.3 

4U 1915-05 

atoll 

655 

1005 

348±11 

3.5 

7.2 

XTE J2123-058 

atoll 

845 

1100 

255±14 

3.5 

9.2 


Low (i^i) and high ( 1 ^ 2 ) QPO frequencies with corresponding peak sep¬ 
aration (Az/) observed in Z and atoll sources. Only one data series is 
given as example, for further detail see van der Klis (2000). Values of 
Ms/R^ and fis are calculated from Eq. (34) for a fixed value a: — 10 as 
an example. Here we assumed M — 10^^ g s“^ and Rs — 10 km. We 
note that larger frequencies may occur at smaller distances, x < 10. 


z„ 



Figure 1. A schematic and idealized side view of an accreting 
neutron star-disk system. The neutron star’s strong gravity 
causes a very high velocity flow toward the magnetosphere. As 
a result, the magnetosphere is pushed inward in the disk plane 
but balloons outward in direction away from the disk plane. 
Some of the plasma may leave the disk and flow along the field 
lines. The magnetic star-disk connection by a dipolar magne¬ 
tosphere as well as the polar accretion flow along the held lines 
are part of a reasonable but simplified picture of magnetized 
star-disk systems in general. The dashed box represents the 
model that we considered in this paper. Zmax > 0 represents 
the surface of the star, Zmin > 0 is on the disk plane, Xmin and 
Xmax are embeded the density gradiant region. The directions 
of the ambient magnetic field and the plasma flows in the box 
are indicated by the verticle arrows, while the direction of the 
gradient in the Alfven speed is shown by the horizontal arrows. 
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r/R^ 


Figure 2. Radial profile of the relativistic Alfven speed va (in 
units of speed of light) in the equatorial plane, see Eq. (30). 
Rs is the radius of the neutron star. Here we assumed p = 0, 
Rs = 10 km, M = Mq, M = 10^^ g s“^ and /r = 10^® G cm®. 



r/R^ 


Figure 3. Frequency of the regular Alfven resonance (vp = 0) 
as function of r, uja = kzVA{R 3 )(r/Rs)~^^'’‘■ The closer to the 
star the higher the frequency. From bottom to top ua(Rs) = 
.2c, .4c, .6c and .8c, respectively. Here Rs = 10 km, and kz = 

•k/Rs- 



Figure 4. Frequency of the regular Alfven resonance {vp = 0) 
as a function of the Alfven wave speed at the stellar radius, 
uja = kzVA{Rs)[r/Ra)~^^'^■ The larger the Alfven speed the 
higher the frequency. From top to bottom r/Rs = 1, 2, 3, and 
4, respectively. Here Rs = 10 km, and kz = tt/R s- 



Figures. Variation of the spectral peak separation Au with 
the lower frequency see Eqs. (31a) and (31c). The result¬ 
ing curves are compared with the observed values for dif¬ 
ferent systems. The solid curves from bottom to top cor¬ 
respond to kzVp{Rs) ~ 650, 750, and 860 Hz. The dashed 
curve corresponds to kzVp{Rs) ~ 710 Hz. Here we assumed 
fJ -26 = .32, Rs = 10 km, Mir = 1, and M — Mq. 
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Figure 6. Variation of the spectral peak separation Av with 
the upper frequency z/ 2 , see Eqs. (31b) and (31c). The re¬ 
sulting curves are compared with the observed values for dif¬ 
ferent systems. The solid curves from bottom to top cor¬ 
respond to kzVp{Rs) ~ 650, 750, and 860 Hz. The dashed 
curve corresponds to kzVp{Rs) ~ 710 Hz. Here we assumed 
M 26 = -32, Rs = 10 km, Mn = 1, and M = Mq. 







